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Abstract. The NA60 experiment has studied low-mass muon pair production in In- In collisions at 158AGeV 
with unprecedented precision. With these results there is hope that the in-medium modifications of the 
vector meson spectral function can be constrained more thoroughly than before. We investigate in particular 
what can be learned about coUisional broadening by a hot and dense medium and what constrains the 
experimental results put on in-medium mass shift scenarios. The data show a clear indication of considerable 
in-medium broadening effects but disfavor mass shift scenarios where the p-meson mass scales with the 
square root of the chiral condensate. Scaling scenarios which predict at finite density a dropping of the 
p-meson mass that is stronger than that of the quark condensate are clearly ruled out since they are also 
accompanied by a sharpening of the spectral function. 

PACS. 25.75.-q Relativistic heavy-ion collisions 
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1 Introduction 

The properties of vector mesons, especially the p-meson 
are expected to change in the a hot and dense nuclear 
medium as created in ultra-relativistic heavy ion collisions. 
The CERES and HELIOS collaborations studied dilepton 
production in central nucleus-nucleus collisions and ob- 
served a strong enhancement of the emission of low-mass 
dileptons compared to scaled proton-nucleus and proton- 
proton collisions 1,2,3,4 . 

It has been argued that these experimental findings 
could be an indication of chiral symmetry restoration. 
The approximate chiral symmetry of Quantum Chromo- 
dynamics (QCD) is spontaneously broken in the vacuum 
but expected to be restored at high temperatures or densi- 
ties. Indeed lattice QCD calculation predict that the prop- 
erties of strongly interacting matter change drastically at 
a temperature about Tc « 170MeV above which chiral 
symmetry is restored and quarks and gluons are decon- 
fined. These calculations also indicate that this transition 
is (at zero net baryon density) a rapid crossover from the 
low-temperature hadron dominated phase into the high- 
temperature plasma phase. For a review, see e.g. 0. 

A gradual restoration of chiral symmetry should also 
lead to modifications in the hadronic sector (at high den- 
sities and temperatures) , since chiral partners become de- 
generate once chiral symmetry is fully restored. Chiral 
symmetry especially dictates a relationship between the 
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vector and axial- vector channels [HIIZI- This is important 
for the dilepton production in a medium since the vector 
correlator dominates the electromagnetic spectral func- 
tion. Therefore it is expected that chiral symmetry restora- 
tion should manifest itself in modifications of the proper- 
ties of the p, Lu and (f> which leave its trace in the measured 
dilepton spectra. 

Brown and Rho proposed that chiral symmetry restora- 
tion should primarily be reflected in a scaling of the masses 
of the vector mesons in a nuclear medium with the quark 
condensate fSnQwlOiall. . This would cause the maximum of 
the yO-meson spectral function to be shifted to lower masses 
in a hot and dense nuclear medium and thus would lead 
to an enhancement of the lepton pairs produced at lower 
invariant masses as compared to the vacuum. 

We also mention that Brown-Rho scaling is not the 
only conjecture that indicates substantial mass shifts. A 
different approach is based on the assumption that the 
lower-mass meson multiplets should be the relevant de- 
grees of freedom near the phase boundary. Pisarski pro- 
posed to study a gauged linear cr-model to gain insight 
into vector meson modifications at finite temperature |12l 
^3 ,14 . He found in a finite temperature mean-field cal- 
culation assuming vector meson dominance that the p- 
meson and ai-meson masses become degenerate at Tc re- 
sulting in rrip — nia^ — 962 MeV. This result is sensitive 
to the vector-meson dominance conjecture, since a situa- 
tion where vector meson dominance does not hold leads 
to slightly decreasing masses close to the phase transition 



Tc, namely nip 



629 MeV. 
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An alternative picture to these scenarios, based on 
broadening of the p-meson without a substantial mass 
shift of the p-meson due to interactions with the surround- 
ing thermal medium was suggested by the study of effec- 
tive Lagrangians in hadronic model calculations (see e.g. 
[TSlllSlllTp . These findings are also supported by calcula- 
tions of the spectral function as inferred from experiment 
by constructing scattering amplitudes for vector mesons 
scattering from pions and nucleons |18lll9j . This substan- 
tial broadening without strong mass shifts of the vector 
mesons can also be in accord with the gradual restoration 
of chiral symmetry if the axial- vector channel exhibits con- 
gruently features. 

The mass resolution of the dilepton spectra of the 
CERES experiment were not high enough to conclusively 
decide which of the scenarios is realized in nature 15 . The 
high precision di-muon data for 158AGeV In- In collisions 
as measured by the NA60 collaboration 20 now provide a 
possibility to put stronger constraints on our understand- 
ing of in-medium effects on vector-mesons. Towards this 
aim we conduct a comparison of different model calcula- 
tions with data. 

Our strategy is as follows: Ideally, one would need to 
construct an evolution scenario for In-In collisions that is 
directly related to hadronic observables in In-In collisions, 
however in practice these are not yet sufficiently well de- 
termined experimentally. This is why we rely on scaling 
arguments to make a sophisticated attempt at describing 
the evolution of matter in 158AGeV In-In collisions based 
on an evolution model for Pb-Pb collisions which has been 
successfully used to describe the plethora of observables 
in Pb-Pb collisions |21| . 

Using this evolution model, we then calculate the re- 
sulting dilepton spectra for different scenarios of in-medium 
modifications of vector mesons and compare with data. We 
are specifically interested in what the NA60 data can tell 
about dropping mass scenarios and broadening effects. 



2 The evolution model 

A detailed description of the evolution model is found in 
j21| . We restrict ourselves to a brief discussion here. The 
scaling from 158 AGeV Pb-Pb to In-In collisions is dis- 
cussed in more detail in 

The evolution model's main assumption is that an equi- 
librated system is created at a short time tq after the on- 
set of the collision. Furthermore, it is assumed that the 
thermal fireball expands isentropically until the mean free 
path of the particles exceeds the dimensions of the systems 
(at a freeze-out time Tf). 

We use the following ansatz for the entropy density: 



distributions 



s{T,r]s,r)^NR{r,T)-H{^s,T) 



(1) 



with T the proper time as measured in a frame co-moving 
with a given volume element and R(r,T), H(r]s,T) two 
functions describing the shape of the distribution and TV a 
normalization factor. They are chosen to be Woods-Saxon 



R{r,T) = 1/ (^1 + exp 
H{r],,T) = 1/ (l + exp 



r - R,{t) 
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to describe the shapes for a given t. 

The ingredients of the model are therefore the skin 
thickness parameters d^s and Ty^s and the parametrizations 
of the expansion of the spatial extensions Rc{t), Hc{t) as a 
function of proper time. For a detailed account how those 
are parameterized see '22'''23j. 

The parameters of the model have been adjusted to 
hadronic transverse mass spectra, HBT correlation mea- 
surements in 158AGeV Pb-Pb and Pb-Au collisions at 
SPS. Since at this point a freeze-out analysis or HBT 
correlation data are not available from NA60 we rely on 
scaling arguments. 

The change in the total entropy production can be 
inferred from the ratio of charged particle rapidity den- 
sities dNch/dr] at 30% peripheral Pb-Au collisions with 
2.1 < 77 < 2.55 measured by GERES PIISIB| and semi- 
central In-In collisions at 77 = 3.8 measured by NA60 j2()j . 

The number of participant baryons and the effective 
initial radius are obtained with nuclear overlap calcula- 
tions. We assume that the stopping power scales approx- 
imately with the number of binary collision per partic- 
ipant, Nfjin/Npart- The shape parameters dws and rj^s 
have no great impact on electromagnetic emission into the 
midrapidity slice, they primarily govern the ratio of sur- 
face to volume emission of hadrons. Thus, we leave them 
unchanged from their value determined in |21j . Under the 
assumption that the physics leading to equilibration is pri- 
marily a function of incident energy we keep the formation 
time r = 1 fm/c as in Pb-Pb collisions. We stress that the 
final results exhibit no great sensitivity to the choice of 
either 770 or tq except in the high M, high pT limit. 

The largest uncertainty is the choice of the decoupling 
temperature Tp. Due to the smaller system size of In- 
In, a higher decoupling temperature can be expected. An 
unambiguous answer could be obtained from HBT corre- 
lations and the transverse mass spectra, for the time be- 
ing we tentatively choose Tp = 130MeV. This translates 
primarily into an uncertainty in the normalization of the 
results. The flow velocity is fixed at Tc by choosing the 
value of transverse flow obtained at T = ISOMeV in the 
Pb-Pb model in the same centrality class. 

There is likewise an uncertainty (which is to some de- 
gree correlated with the decoupling temperature) since 
flow determines the slope of hadronic (and to a lesser de- 
gree electromagnetic) pr spectra. 

The equation of state in the hadronic phase and the 
off-equilibrium parameters ^tt^^ik are inferred from sta- 
tistical model calculations as described in |21II24| . The re- 
sulting fireball evolution is characterized by a peak tem- 
perature of about 250 MeV, a lifetime of ^ 7.5 fm/c and 
a top transverse flow velocity of 0.5c at decoupling. 
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3 The hadronic spectral function 

Different effective liadronic models and techniques for the 
calculations of the properties of hadronic matter near and 
below the phase boundaries have been suggested, for re- 
views of the intense theoretical activities, see e.g. [TKHT?)! 
[T7j . Most of those calculations predict substantial broad- 
ening in matter with comparably small shifts of the in- 
medium masses. However, this might be different if one 
follows the conjecture by Brown and Rho that in-medium 
vector- meson masses decrease considerably due to the grad- 
ual restoration of chiral symmetry and a scaling of the 
masses with the chiral condensate in the vicinity of the 
phase transition |51i^lOllll| . We therefore discuss this is- 
sue separately from hadronic model calculations. 

3.1 Hadronic model calculations 

We briefly discuss three different hadronic model calcu- 
lations which predict broadening effects but small mass 
shifts as generic examples. The first model is based on 
j25ll26| where the electromagnetic current-current correla- 
tor has been computed in the so called improved vector 
meson dominance model combined with the chiral dynam- 
ics of pions and kaons. We refer to this model as 'chiral 
model' further on. For the evaluation of finite baryon den- 
sity effects which are relevant at SPS conditions the re- 
sults from Klingl, Kaiser, and Weise j27j are used. Ther- 
mal broadening of the p, and w were calculated by 
Schneider and Weise as modifications on top of the finite 
baryon density effects using perturbative methods 26 , 25 . 
In these calculations it is assumed that the temperature 
and density dependences of the vector meson self-energy 
factorize. This amounts to neglecting contributions from 
pion-nucleon scattering where the pion comes from the 
heat baths. Formally speaking matrix element contribu- 
tions such as {-kNXT i^(x)i>^{x)\'KN) are neglected. 

A second approach ^22, (referred to as 'Walecka model' 
in the following) is based on a different idea: The method 
of thermofield dynamics is used to calculate the state 
with minimum thermodynamic potential at finite temper- 
ature and density. The temperature and density depen- 
dent baryon and sigma masses are calculated selfconsis- 
tently. The medium modification to the masses of the lo- 
and p-mesons in hot nuclear matter including the quan- 
tum correction effects are then calculated in the relativis- 
tic random phase approximation. The decay widths for 
the mesons are calculated from the imaginary part of the 
self energy using the Cutkosky rule. 

We also discuss a third model calculation where we 
solved truncated Schwinger-Dyson equations incorporat- 
ing a self-consistent resummation of the tt — p interaction 
in order to determine the thermal broadening of the p- 
meson, see [211 and especially note [211. We refer to this 
model as the tt — p model. This ^-functional approach 
self-consistently takes into account the finite in-medium 
damping width of the pion in the thermal heat bath. Ef- 
fects of baryons are not yet included in the self-consistent 
approximation scheme. The solution of the ^-functional 



approach includes the full momentum dependence of the 
three-dimensional longitudinal and transverse components 
of the p-meson spectral function which is especially im- 
portant for a comparison with the NA60 data measured 
at integrated and high px, see section II D. 

3.2 Mass scaling model 

The universal scaling of the masses of hadronic particles, 
especially for the p-meson with the chiral condensate were 
proposed by Brown and Rho Given this conjecture, 
the gradual restoration of chiral symmetry would cause 
the maximum of the spectral function of the p-meson to 
be shifted to lower masses in heavy-ion collisions. This 
would lead to an enhancement of lepton pair production 
at lower invariant masses as compared to the vacuum. 
Calculations in the hidden-local-symmetry model at finite 
temperature indeed support the notion of dropping masses 
31,32 . But an explicit calculation of spectral functions in 
such an approach, including finite temperature and baryon 
density as well as momentum dependence of the spectral 
function is not yet available. 

Therefore we take a different approach here and re- 
strict ourselves to the discussion of simplified schematic 
scaling models. 

We construct a simple mass scaling model which allows 
for the investigation of a dropping of the p-mass scaled 
with the quark condensate and takes in-medium changes 
of the width of the p-meson in a nuclear medium into 
account. 

This model has three main assumptions: i) the relative 
change of the in-medium mass of the p-meson is propor- 
tional to the relative drop of the quark condensate to some 
power a in the medium: 



(ii) the width of the p-meson can change in the medium, 
and (iii) the validity of vector meson dominance is as- 
sumed. For simplicity we assume the rho-meson to be at 
rest and employ a Breit-Wigner parametrization for the p- 
meson spectral function ■33"'34! and note that this ansatz 
may not account for peculiar features, e.g. two separate 
pole structures, that could appear in certain classes of 
hadronic models. 

Taking into account the change of the width of the 
p-meson is necessary since predictions for the in-medium 
mass and in-medium width of the p-meson in a nuclear 
medium cannot be independent but have to be strongly 
correlated in order to fulfill the QCD sum rule approach 

For simplicity we restricted our analysis to QCD sum 
rules at finite baryon density only j34| and did not include 
temperature effects on the condensate/spectral function. 
Although within the (T, pb) trajectory probed by the evo- 
lution in the In-In system baryon density seems to be 
the leading influence on the condensate, our model should 
only be taken as a lower limit of the potential downward 
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mass scaling effect, the actual mass shift realized in certain 
space-time regions of the collision may be stronger. 

In |34| it was shown that the lower and upper bounds of 
the in-medium width of the p-meson at rest that are com- 
patible with the QCD sum rules depend on how rapidly 
the mass decreases in comparison to the change of the 
quark condensate. The bounds for the in-medium width 
increase with density if the relative change of the quark 
condensate is stronger than the relative decrease in mass 

Assuming a specific scaling law one can extract the 
upper and lower bounds of the in-medium width of the 
/9-meson at a given density. Those bounds are sensitive to 
the parameter sets for the QCD sum analysis Dif- 
ferences between the bounds for different parameter sets 
show intrinsic uncertainties of the approach. Using these 
boundaries for the in-medium width one can study the 
question if the NA60 data are compatible with a specific 
scaling law of the in-medium p-meson mass. 

To be specific, we consider here the case a=l/2 in [S^ 
(see Fig. 1 therein). We analyzed the upper and lower 
bounds of compatible width in this case for three different 
parameter sets of the QCD sum rule approach '35 , 27' 36 . 
The case a = 1/2 has some resemblance to the Brown- Rho 
scaling law |10llll| if one identifies the parametric mass of 
Brown- Rho scaling with the in-medium p-meson mass |37l 
I38j . This could be a reasonable approximation up to about 
normal nuclear matter density, ps — po, since the dense 
loop term contribution is expected to be small jSZIEH!- In 
spite of these obvious similarities between our mass scal- 
ing scenario and the Brown- Rho scaling law one should be 
cautious and consider the arguments by Brown and Rho 
|37II38| why an identification of Brown-Rho scaling with 
simplified mass scaling models is problematic. They pro- 
pose to address the scaling issue (especially also the iden- 
tification of parametric and in-medium p-meson mass) in 
the framework of the hidden local symmetry with vector 
manifestation pT^31ll32j . Although the issue of a possible 
violation of vector meson dominance has been stressed re- 
cently, this question cannot be expected to be of primarily 
importance for an understanding of the lineshape of the 
spectra but would effect the overall normalization. 

4 Results 

We first present a comparison with the data using the 
spectral function for the tt — p systems for all px- The cal- 
culation contains a vacuum p contribution using a Cooper- 
Frye prescription (for details see section II E of [23 )• Most 
of the vacuum p-contribution is generated during the fi- 
nal breakup of the fireball j22ll23| and thus the underlying 
flow is strong and shuffles strength to higher pT- This is 
the reason that the vacuum p-contribution is almost neg- 
ligible for the pt < 0.5 GeV spectra, but is important for 
integrated pt (and high px > IGeV). The relative con- 
tribution of the vacuum p and the in-medium p appears 
to be approximately correct. While the strength of the 
vacuum p parametrically scales like the freeze-out volume 
times the p-meson's life-time, the in-medium contribution 



is sensitive to the lifetime of the fireball. Thus the good 
description in the mass region M « 0.6 — 0.9GeV is due to 
contributions of the vacuum and in-medium p and there- 
fore a consistency check of the fireball evolution. We note 
that the pr > IGeV spectra are in good agreement with 
the NA60 data (comparison not shown). 

The good agreement of the <?-derivable approach cal- 
culation with the data including the vacuum p is to some 
degree surprising and unexpected, as the model contains 
the self-consistent broadening due to the tt — p interac- 
tion ^29.30 but is completely insensitive to the presence 
of other hadrons, especially baryons (which, based on per- 
turbation theory, are expected to play a major role in the 
broadening of the p). One may take this as an indication 
that the generic shape of the spectral function and its 3- 
momentum dependence is right, not as a proof that this is 
the main mechanism capable of producing such a spectral 
function or that the dynamics of baryons is unimportant. 
However, these results clearly indicates that it is neces- 
sary to extend studies on non-perturbative mechanisms 
for broadening in the future. 

The other model calculations (chiral model, Walecka 
model, mass scaling scenarios) rely on the assumption of 
a p-meson at rest relative to the thermal heat bath q = 
0, hence they can only be compared to the data in the 
Pt < 0.5GeV cut, see This is somewhat unfortunate 
since the acceptance of the NA60's experimental setup is 
relatively small in the M < 0.6 MeV, pt < 0.5GeV region, 
nonetheless interesting conclusions can be drawn. 

The chiral model and the Walecka model in Fig. ^ 
give a fair description of the data in the low px region. 
Although we note that the chiral model calculation also 
contains an in-medium and lu contribution (as these have 
different mass/width in medium than in vacuum such a 
contribution would not be removed by the experimental 
subtraction of the cocktail 4>, uj) and apparently the strong 
in-medium (j) broadening is not in good agreement with the 
low-pT data. 

We present the results of the calculation with the mass 
scaling scenarios as described in 13.21 in Fig. ^ after ac- 
ceptance folding and maintaining the normalization for 
the two different parameter sets PI [21] and P2 |27j used 
in the QCD sum rule approach "33",''34' . Since the spec- 
tral function is determined for vanishing three momentum, 
only comparison with the low px binned data is meaning- 
ful. Given this comparison one realizes that such a p-mass 
scaling behavior is disfavored in the face of the experi- 
mental data. Although part of the observed enhancement 
below the mass of the free p-meson is accounted for, the 
yield in the peak region is underpredicted. As discussed 
above the vacuum p-contribution is not strong enough to 
account for this difference. What also can be stated is that 
scaling scenarios which would predict a dropping mass 
of the p-meson with the quark condensate that is much 
stronger than in the a = 1/2 case, i.e. a ^ 1, are in dis- 
agreement with the data. This has two reasons, namely, (i) 
they predict a stronger decrease of the p-meson mass with 
density leading to a stronger shift of the "peak" region 
down to small invariant masses of the di-muon spectra. 
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Fig. 1. Left panel: Comparison of the NA60 dimuon data |2()) with calculations within a self-consistent functional approach 
for the n — p interaction '29, with the all pr-data. Right panehComparison of the NA60 dimuon data |2U| within a (^-functional 
approach for the n — p interaction (RR) 1291 . a chiral model ('RSW/KKW') j25) and a calculation in the Walecka model (MR) 
|28| . Model calculations for the chiral model and Walecka model have been normalized to the data in the mass interval Af < 0.9 
GeV. 




Fig. 2. Comparison of the NA60 dimuon data '20' for pr < 0.5 GeV with schematic model calculations of Brown-Rho scaling 
in combination with the broadening compatible with QCD sum rules as described in |34| . Left panel: lower limits of the width 
for two parameter sets PI to P2 describing the density dependence of the condensate. Right panel: upper limits of the width 
as determined by sum rules. The 3-momentum dependence of the spectral function has been neglected. All model calculations 
have been normalized to the data in the mass interval M < 0.9 GeV. 



and (ii) a stronger drop of the mass demands a much less 
rising or for a > 1 even decreasing in-medium width of the 
p- meson with density |34| . both in disagreement with the 
NA60 dimuon data. We note that Rapp and van Hees have 
studied a mass scaling model PHI recently, where they use 
a mass parameterization of the type 

m; - mp(l - 0.15pb/po)[1 - {T/T,)^' (4) 

and employed vector meson dominance. They supplemented 
the in-medium spectral function by thermal broadening 
and obtained qualitatively similar results for mass-scaling 
scenarios, namely somewhat an underprediction of the 
yield in the peak region due to shift to lower masses. We 



mention that Greiner and Schenke recently argued that 
non-equilibrium effects can lead to significant different re- 
sults in comparison to equilibrium calculations due to the 
importance of memory effects 40- They find that this 
might have significance for the normalization and line- 
shape of dilepton spectra as obtained in mass scaling sce- 
narios |41| . 

5 Conclusion 

In these proceedings we have shown how a dynamical evo- 
lution model of 158 AGeV In-In collisions can be con- 
structed in order to calculate the measured di-muon yield 
by applying scaling arguments to an evolution model that 
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is able to describe the plethora of observables in Pb-Pb 
collisions |21j . 

Furthermore we discussed how the low-mass NA60 di- 
muon measurements can be used to disentangle different 
scenarios: strong broadening of the spectral function with 
a small shift of the p-meson mass in medium as inferred 
by hadronic-many body calculations vs. considerable mass 
shift of the p-meson. We have shown that taking into ac- 
count only a mass shift without considering changes of the 
in-medium width in a nuclear medium violates QCD sum 
rule constraints ,34, ■ We therefore presented a schematic 
mass scaling model in order to test mass scaling assump- 
tions in a dense medium taking the correlation between 
mass shifts and finite width effects into account. To ac- 
count for the intrinsic uncertainties we used two different 
parameterization in the QCD sum rule approach. It was 
shown that the evolution model together with hadronic 
many body calculations which include considerable broad- 
ening can account for the data whereas the same evolution 
model assuming a mass scaling scenario underpredicts the 
yield in the " peak" region of the data and indicates a shift 
of the maximum of the di-muon spectrum down to smaller 
invariant masses. 

The NA60 data therefore do not favor downward mass 
scaling and clearly indicate that the dominant in-medium 
effect in the medium created in 158AGeV In-In collisions 
is broadening of the p-spectral function. 
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